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Exact electronic properties of a system of four Coulomb- 
interacting two-dimensional electrons in a parabolic confine- 
ment are reported. We show that degenerate ground states 
of this system are characterized by qualitatively different in- 
ternal electron-electron correlations, and that the formation 
of Wigner molecule in the strong-interaction regime is going 
on in essentially different ways in these ground states. 

PACS numbers: 73. 21. La, 71.10.-w, 73.20.Qt 



Quantum dots - artificial electron systems real- 
izable in modern semiconductor structures, - are ideal 
physical objects for studying effects of electron-electron 
correlations. They contain a few (N) two-dimensional 
(2D) interacting electrons moving in the plane z = 
in a lateral confinement potential V(x,y). The poten- 
tial V can often be modeled by a harmonic oscillator 
form V = K(x 2 + y 2 )/2, with a characteristic frequency 
luq = y/K/m* and a characteristic length Iq = y/K/ m*u>o 
(here m* is the effective mass of electrons) . The relative 
strength of interaction and kinetic effects can be charac- 
terized by the dimensionless parameter 



A = 



7«B 



huj n 3/2 



(i) 



determined by the ratio of quantum-mechanical and in- 
teraction length (energy) scales {as is the effective Bohr 
radius of the host semiconductor) . The curvature of the 
confinement K, and hence the interaction parameter A, 
can be experimentally varied, so that quantum dots are 
often referred to as artificial atoms with tunable physical 
properties. Theoretical study of physical properties of 
such "atoms" as a function of the Coulomb-interaction 
parameter A may give a valuable information about the 
role of electron-electron correlations in interacting quan- 
tum systems. The physics of interactions becomes espe- 
cially interesting in zero external magnetic field B, when 
electronic spins are not polarized and are active play- 
ers in the game. Detailed theoretical study of physical 
properties of quantum-dot atoms, including the Fermi- 
liquid - Wigner-molecule crossover in the ground state 
with growing strength of intra-dot Coulomb interaction, 
is one of the challenging problems in this field, which was 
attracting increasing interest in recent years [PHf5| (for 
more references see [[uj and a recent review jig] ). 
The Hamiltonian of an ./V-electron dot 



2m* 



1 N 



(2) 



commutes with the operators of the total angular mo- 
mentum L = L z , total spin S, and its projection on 
some axis. Solutions 



E n i = hu (2n 



1) 
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of the single-particle problem of quantum-dot hydrogen 
(N = 1) M] are characterized by the radial n and 
asimutal I (angular momentum) quantum numbers. The 
ground state of quantum-dot helium [N = 2) [|18j is char- 
acterized by the total angular momentum L = and 
the total spin S = 0. In the limit of weak interaction 
A <C I two electrons with opposite spins occupy the 
single-particle state (n,l) — (0,0), forming the config- 
uration (0,0, tl). The state (L,S) = (0,0) remains the 
ground state at any A, also in the regime of strong inter- 
action. The first excited state of quantum-dot helium is 
the state (L, S) = (1, 1). 

In quantum-dot lithium (N = 3) |Q ] the A-dependence 
of the ground-state properties is more complicated. In 
the regime of weak interaction A <C 1 electrons form the 
configuration [(0, 0, Tl), (0, 1, T)] with (L, S) = (1,1/2). 
The state (1, 1/2) remains the ground state up to A = 
A c w 4.343, where a transition to another ground state 
with (L,S) = (0,3/2) occurs. All physical properties of 
the dot sharply change at the transition point. At A > A c 
electron distribution in the dot has the form of an angle- 
averaged equilateral triangle (the Wigner molecule). 

The goal of this Letter is to present exact results of so- 
lution of the four-electron quantum-dot beryllium prob- 
lem at B = 0. Some results for the energy of states at 
a few points on the A-axis, obtained by different vari- 
ational methods, were reported in the literature (e.g. 
P,[7|JT2]|). We complete and essentially improve these re- 
sults, presenting them for the ground and a number of 
excited states in a broad range and as a function of A. 
A principally new finding of this work results from a de- 
tailed study of ground-state properties of quantum-dot 
beryllium. The ground state of this atom is degener- 
ate with respect to the total spin projection and we 
show that electron-electron correlation functions, as well 
as their spin density, are substantially different in dif- 
ferent ground states. As a result, formation of Wigner 
molecules with the growing interaction strength is going 
on in qualitatively different ways, dependent on the value 
of the total spin projection St (0 or ±1). This result is 
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valid not only in the considered parabolic quantum-dot 
model. It has a general character and should be the case 
in any strongly-interacting many-electron system (e.g. in 
natural atoms) with degenerate ground states. 

The Schrodinger equation with the Hamiltonian (j^) 
and N = 4 is solved by the exact-diagonalization method 
described in Ref. [jlij. The many-body wave function 
is expanded in a complete set <I> M of many-body eigen- 
states of the Hamiltonian Hq of non-interacting parti- 
cles, H^ u = E u ^ u , with the number of basis states 
N m bs restricted by the condition E u < E t h, where E t n 
is a threshold value. The larger the threshold energy 
E t h, the higher the accuracy of the calculated energy at 
a given value of A, and the broader the range of A's, in 
which results have a given accuracy. We performed cal- 
culations until a sufficient (very high) accuracy has been 
achieved in the whole (studied) range of A (A < 10). For 
the ground state energy our accuracy is about two or- 
ders of magnitude higher than in corresponding Quantum 
Monte Carlo (QMC) calculations; at A = 10 it comprises 
about 6 x 10~ 4 % (see Table II below). 

Figure [l]a shows the ground-state energy of a four- 
electron parabolic quantum dot as a function of the in- 
teraction parameter A, together with the classical contri- 
bution 

E cl = 6[(1 + VS)/4] 2 / 3 (e 2 /l cl ), l cl = (e 2 /m*^)i/3 ) 

(4) 

(the energy of four Coulomb-interacting point-like parti- 
cles in a parabolic potential fl9l , |20|| ) , and with the curve 
E c i + 4?iwo (the classical energy plus htuo per particle). 
The ground state is characterized by quantum numbers 
(L, S) = (0, 1) and is three-fold degenerate. In the weak- 
interaction regime the corresponding wave functions are 
*i, (* 2 + *3)/V2 and * 4 (S c = +1,0, and -1, respec- 
tively), where 

*i = [(o,o,n),(o,-i,T),(o,i,T)L 
*2 = [(o,o,n),(o,-i,T), (o,U)], 
*3 = [(o,o,u),(o,-i,|),(o,i,T)], 
* 4 = [(0,0, u), (o,-i, 4), (o,i, l)]. (5) 

Figure 0b shows the energies of a few excited states (the 
ground-state energy is subtracted). At small A the first 
and the second excited states are (L,S) = (2,0) and 
(L, S) = (0, 0) respectively. The state (L, S) = (2, 0) is 
two-fold degenerate and corresponds, in the limit A —> 0, 
to the configurations [(0, 0, U), (0, -1, T4)] and [(0,0, U 
), (0,l,t|)]- The state (L,S) — (0,0) is non-degenerate 
and corresponds to the configuration (ty 2 — ^ , 3)/v / 2 at 
A — ► 0. These and some other excited states are listed 
in Table | in the order of increase of their energies at 
A< 1. At larger A's the order of excited levels changes, 
as seen from Figure |l|b. The level crossings are the case at 
A w 3.40 [the states (L,S) = (0,0) and (2,0)], A « 4.42 



[the states (1,1) and (2,0)], and A « 2.37 [the states 
(2,2) and (3,1)]. The state (L,S) = (0,1) remains the 
ground state at all (studied) values of A (< 10). 

Table || exhibits numerical data, at a few A points, 
for the ground state (L, S) = (0, 1), and for one of the 
excited (fully spin-polarized) state (L, S) — (2,2). Here 
we also compare our results with QMC calculations from 
Ref. 0- Notice that, although the QMC results are, in 
general, in a good agreement with exact ones for the fully 
spin-polarized (excited) state, they have an essential sys- 
tematic error for the ground state (L,S) = (0, 1): for in- 
stance, at A = 10 the absolute error of the QMC method 
(about 0.07?kJo) is substantially larger than the energy 
difference between the states (0, 1) and (0, 0) (estimated 
as < 0.013?iwo at A = 10). 

Now consider physical properties of quantum-dot 
beryllium in the ground states (L,S,S^) — (0,1,0) and 
(L, S, S{) = (0, 1, +1). These states are essentially differ- 
ent. In the ground state with = two electrons are 
polarized up and two electrons are polarized down. They 
form a symmetric square structure, with two up- and two 
down-polarized electrons occupying opposite corners of a 
square, see Figures || and || showing the density and the 
pair-correlation functions of electrons in this state. The 
relative electronic distribution holds its square form at 
all A's, also in the regime of small interactions A <C 1. In 
the ground state with = +1 three electrons are polar- 
ized up and one electron is polarized down. In the weak- 
and moderate-interaction regimes (A < 2) they form a 
triangular structure with one (down) electron occupy- 
ing the center of the dot, and three (up) electrons rotat- 
ing around, see Figures | and |. When A increases, the 
growing Coulomb interaction pushes the down-electron 
out from the center, and the triangular configuration is 
smoothly transformed to a (classical-type) square-form 
structure. In Figure || one sees how this process begins 
at A ~ 3 (the third row of plots) and how it is completed 
at A ~ 4.5 (the last row). Thus, correlations between 
electrons with different spins are essentially different in 
ground states with — and = ±1, and the for- 
mation of Wigner molecules is going on in qualitatively 
different ways in these states. 

The coexistence of degenerate ground states with qual- 
itatively different electron-electron correlations is obvi- 
ously independent of the 2D nature of the electrons and 
of the parabolic form of the confinement potential in the 
present model. The same situation should be also the 
case, for example, in a natural carbon atom, which has 
S = 1 in the ground state. In a hypothetic case of a quan- 
tum system with S = 2 and more particles one should 
also expect that degenerate ground states with = 0, 
±1 and ±2 are characterized, in the strong-interaction 
regime, by qualitatively different inter-electron correla- 
tions. 

To conclude, we have studied energy spectra, charge 
and spin densities, and electron-electron correlations in 
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quantum-dot beryllium - a system of four strongly inter- 
acting electrons in a harmonic oscillator potential. It was 
shown that this artificial atom can exist in two essentially 
different ground states, characterized by qualitatively dif- 
ferent electron-electron correlation functions. 

The work was supported by the Deutsche Forschungs- 
gemeinschaft through SFB 484. 
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FIG. 1. (a) The energy of the ground state (L, S) = (0, 1) 
of quantum-dot beryllium as a function of the interaction pa- 
rameter A, together with the clasical contribution E c i and the 
energy E c i + 4Tiljo- (b) The energies of a few excited states 
E(ls) — ^(oi) (the ground state energy subtracted) as a func- 
tion of A. Calculations were performed with E t h = 22hui , 
N m bs = 14453. The energy unit on both plots is huio- 
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FIG. 2. Electron density 7rZon CT (r) of spin-up and 
spin-down polarized electrons [ni(r) = ni(r)] as a function 
of radial coordinate r/lo in the ground state (L,S) = (0, 1) 
with S( = at a few values of the interaction parameter A. 
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FIG. 4. The same as in Figure g, but for the ground state 
with S( = +1. The total (charge) density n-f(r) +nj,(r) is the 
same as in Figure |2| 








FIG. 3. Pair-correlation functions P^(r,r') = P^(r,r') 
(left column) and Pn(r, r') = Pf|(r, r') (right column) in the 
ground state (L, S) = (0, 1) with — 0, as a function of 
r at r' = (0, — R c {) (the length unit is Iq). The interaction 
parameter A assumes the values A = 0.2, 2, and 8, from up to 
down (R c i/lo = 0.57, 1.24, and 1.97, respectively). 
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FIG. 5. Pair-correlation functions Pff(r, r') (left col- 
umn), and P||(r, r') (right column) in the ground state 
(L,S) — (0,1) with S( = +1, as a function of r at 
r' = (0, —Rci) (the length unit is lo, the second subscript 
corresponds to the electron fixed at the r' point). The in- 
teraction parameter A assumes the values A = 0.2, 2, 3, and 
4.5, from up to down [R c i/l = 0.57, 1.24, 1.42, and 1.62, 
respectively) . 

TABLE I. The ground and a few lowest excited states of 
quantum-dot beryllium: the total angular momentum and 
spin (L, S), energy in the weak-interaction limit .Els (A — » 0) 
(in units hu)o), and degeneracy g of the states. The states are 
shown in the order of increasing of their energy at A <C 1. 



No 


(L,S) 




9 


1 


(0,1) 


6 


3 


2 


(2,0) 


6 


2 


3 


(0,0) 


6 


1 


4 


(1,1) 


7 


6 


5 


(3,1) 


7 


6 


6 


(1,0) 


7 


2 


7 


(1,1) 


7 


6 


8 


(1,1) 


7 


6 


9 


(1,0) & (3,0) 


7 


4 


10 


(1,0) 


7 


2 


11 


(2,2) 


8 


10 


12 


(0,2) 


8 


5 



TABLE II. Energies of the ground state (L, S) = (0, 1) and 
of the fully spin-polarized state (L, S) = (2, 2) in quantum-dot 
beryllium from exact diagonalization (this work) and Quan- 
tum Monte Carlo calculations Q . Exact data in the first five 
lines were obtained with Eth = 22Tlujo, which corresponds to 
N m bs = 14453 for the ground state and N m b 3 = 3248 for the 
spin polarized state. The data in the three last lines were ob- 
tained with Eth = 24hLu (Nmbs = 24348) for the state (0, 1) 
and with Eth = 2Qhujo, Nmbs = 8721 for the state (2, 2). 



A 


(0,1) 


(o,i) a 


(2,2) 


(2,2) a 


2 


13.6180 


13.78(6) 


14.2535 


14.30(5) 


4 


19.0323 


19.15(4) 


19.3565 


19.42(1) 


6 


23.5958 


23.62(2) 


23.8025 


23.790(12) 


8 


27.6696 


27.72(1) 


27.8203 


27.823(11) 


10 


31.4122 


31.48(2) 


31.5326 


31.538(12) 


10 


31.4120 


31.48(2) 


31.5323 


31.538(12) 


15 


39.8163 




39.8970 




20 


47.4002 




47.4013 
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